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We study cosmological consequences of non-constant brane world moduli in five dimensional 
brane world models with bulk scalars and two boundary branes. We focus on the case where 
the brane tension is an exponential function of the bulk scalar field, Ub oc exp(Q(/!>). In the limit 
a ^ 0, the model reduces to the two-brane model of Randall-Sundrum, whereas larger values of 
Q allow for a less warped bulk geometry. Using the moduli space approximation we derive the 
four-dimensional low-energy effective action from a supergravity-inspired five-dimensional theory. 
For arbitrary values of a, the resulting theory has the form of a bi-scalar-tensor theory. We 
, show that, in order to be consistent with local gravitational observations, a has to be small (less 

than lO"'^) and the separation of the branes must be large. We study the cosmological evolution 
of the interbrane distance and the bulk scalar field for different matter contents on each branes. 
Our findings indicate that attractor solutions exist which drive the moduli fields towards values 
consistent with observations. The efficiency of the attractor mechanism crucially depends on the 
matter content on each branes. In the five-dimensional description, the attractors correspond to 
the motion of the negative tension brane towards a bulk singularity, which signals the eventual 
' breakdown of the four-dimensional description and the necessity of a better understanding of the 

^ , bulk singularity. 
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I. INTRODUCTION 



One of the most exciting ideas which has originated from theories of particle physics beyond the Standard Model is 
that our universe could be a three-dimensional object (called a brane) embedded in a higher-dimensional space (the 
bulk). Special attention has been paid to models in which the bulk is five-dimensional. For these models, cosmological 
consequences have been worked out in particular for the cases of heterotic M-theory — a five-dimensional theory 
(~| i with two branes, which are the boundaries of the bulk spacetime (see e.g. [1-7]) — and the Randall-Sundrum model 
ILJ ■ [8,9], in both the version with one brane and two branes (see e.g. [10-21]). In the case of the latter model, the bulk 
spacetime is highly warped (in fact, the bulk is an Anti-de Sitter spacetime), whereas in heterotic M-theory the bulk 
is "slightly" warped. Another essential ingredient in heterotic M-theory is the existence of a bulk scalar field, whereas 
^ in the original Randall-Sundrum model there is no such bulk scalar. It was later introduced in order to stabilize the 
inter-brane distance [22,23]. 

It was shown in several papers that the resulting four-dimensional theories at low energies in both models have 
much in common with (multi)scalar-tensor theories, where the interbrane distance and the bulk scalar field degree of 
freedom play the role of scalar degrees of freedom in the gravitational sector of the effective four-dimensional theory 
[24-32]. From the cosmological point of view it is interesting to investigate the consequences of these moduli. Although 
it is usually assumed that these moduli fields are stabilized by some unknown mechanism in the early universe [33-36] , 
it is imaginable that they are not stabilized by some potential. Instead, it might be that there is a cosmological 
attractor mechanism at work, similar to the one found in Brans-Dicke type theories, which drives the theory towards 
four-dimensional General Relativity during the cosmological evolution by generating very small matter couplings of 
the moduli [37-39]. If this is the case, there could be cosmological consequences of the time-evolution of the moduli, 
such as varying constants, time-varying masses, etc., leading to observational consequences, such as in the cosmic 
microwave background radiation, large scale structures, time-varying constants. Equivalence Principle violations, and 
so on. Indeed, there is another motivation for our work: namely, the claim of a time-varying fine-structure constant 
made in [40] and [41]. If this claim is confirmed by future observations, this would suggest that at least some moduli 
fields are not stabilized but slowly varying. 

In this paper we investigate the cosmological evolution of moduli fields arising in brane world models and address 
the question of whether there is a cosmological attractor for the brane world moduli. We also analyse the conditions 
under which conditions such an attractor would be efficient enough for the theory to agree with local experiments in 
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the laboratory and/or strong field limits [42]. We study these issues in a broad class of models encapsulating features 
from both the Randall-Sundrum models and heterotic M-theory, i.e. a warped background with a bulk scalar field. 

The paper is organized as follows: in the next section we discuss the five dimensional theory and derive in detail 
the effective low-energy action from the moduli space approximation, giving also the resulting field equations in the 
Einstein frame. In Section 3 we clarify the conditions under which the theory would predict time-varying coupling 
constants, such as the fine-structure constant, and investigate the constraints imposed by strong field limits. Cos- 
mological considerations and solutions to the field equations are presented in Section 4, and we conclude in Section 
5. 



II. MODULI SPACE APPROXIMATION 



In this section we derive the four-dimensional low-energy effective action using the moduli space approximation. 
For the gravitational sector, our derivation follows closely that in [43]. We begin with the description of the static 
configuration, which was derived from supergravity [44] . 



A. The static configuration 



The bulk action consists of two terms which describe gravity and the bulk scalar field dynamics: 



(1) 



Further, our setup contains two branes. One of these branes has a positive tension, the other brane has a negative 
tension. They are described by the action 

Shra.nel = J d^X^^^^UsSizi), (2) 

J d^X^^^UB5{Z2). (3) 



*^brane2 — ~t~ ^ o 
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In these expressions, Z\ and Z2 are the (arbitrary) positions of the two branes, i/s is the superpotential; U, the bulk 
potential energy of the scalar field, is given by 



We will also include the Gibbons-Hawking boundary term for each brane, which have the form 



= \ I d'^Xy/^K, (5) 



where K is the extrinsic curvature of the individual branes. We impose a Z2-symmetry at the position of each brane. 
The solution of the system above can be derived from BPS-like equations of the form 

a' Ub ., dUs 

— = — T' ^ = ITT' 
a 4 o^) 

where ' = d/dz for a metric of the form 

ds^ = dz'^ + a^{z)r]^,„dx^'dx''. (7) 
We will particularly focus on the case where the superpotential is an exponential function: 

Ub = 4A;e"'^. (8) 
The values a = l/VS, -1/Vl2 were obtained in a theory with supergravity in singular spaces [44] . The solutions read 
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a{z) = (1 - Aka^z)^, 



(9) 



while the scalar field solution is 

tp = --hiil-Aka^z) . (10) 
a 

In the a — > we retrieve the AdS profile 

a{z) = e-'=^ (11) 

Notice that in that case the scalar field decouples altogether. Also, notice that there is a singular point in the bulk 
at = l/Aka^, for which the scale factor vanishes. This will be important for the discussion on the cosmological 
evolution in Section 4. 

In the following we will discuss the moduli space approximation. Wc will put matter on the branes as well as 
supersymmetry breaking corrections to the brane potentials. Two of the moduli of the system are the brane positions. 
That is, in the solution above the brane positions are arbitrary. In the moduli space approximation, these moduli are 
assumed to be space-time dependent. We denote the position of brane 1 with zi = 4'{x'^) and the position of brane 
2 with Z2 = (j{x^). We consider the case where the evolution of the brane is slow. This means that in constructing 
the effective four-dimensional theory we neglect terms like {d(p)^. 

In addition to the brane positions, we need to include the graviton zero mode, which can be done by replacing 77^,^ 
with a space-time dependent tensor (7^^1/(2;^). Thus, we have two scalar degrees of freedom, namely the positions of 
the two branes which we will denote with ^{x^) and (t{x^), and the graviton zero mode g^i, (see [45] and [46]). As 
we will see below, in the coordinate system in which the branes move, the kinetic terms of the moduli come from the 
boundary terms alone. 

It would also be possible to consider another coordinate system, in which the branes are fixed. Then, the moduli 
are a space-time dependent part which can be added to the bulk scalar field solution (10) and the 55-component of 
the metric tensor becomes a four-dimensional effective field. In addition, we have the graviton zero-mode. Thus, the 
number of moduli fields is independent of the coordinate system. 

Note that the moduli space approximation is only a good approximation if the time-variation of the moduli fields 
is small. This should be the case for late-time cosmology well after nucleosynthesis, which we are interested in (in 
[47] and [48] the moduli space approximation was also used in the context of brane worlds). 



B. Moduli space approximation: the gravitational sector 

Let us first consider the bulk action. Replacing 7]^^ with g^^{x'^) in (7) we have for the Ricci scalar R = R^^^/a^+R, 

where R is the Ricci scalar of the backgroimd (7). We explicitly use the background solution (9) and (10), so that R 
will not contribute to the low-energy effective action. Also, in this coordinate system, where the branes move, there 
is no contribution from the part of the bulk scalar field. Collecting everything we therefore have 

5buik = ^/ dzd^xa^y/^^n^^^ = J d^xy/^f{^,a)n^'^\ (12) 

with 

f{^,a) = ^ r dza^z). (13) 

'^5 J<l> 

We remind the reader that a{z) is given by (9). 

We now turn to the boundary terms. It is clear that the integrals (2) and (3) do not contribute to the effective 
action for the same reason that R does not contribute. Let us therefore turn our attention to the Gibbons-Hawking 
boundary terms. 

First, it is not difficult to construct the normal vectors to the brane: 

+ {d(j)y/a^ 

Then the induced metric on each brane is given by 
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ind.l 



Thus, 
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2a2(</,) 
1 



So the Gibbons-Hawking boundary terms take the form 



4I 



j4 4 

a xa 



j4 4 

a xa \ 
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The trace of the extrinsic curvature tensor can be calculated from 

1 



Neglecting higher order terms this gives 



K = 



K 



'-9h 



;d from 



K, 
K. 
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(15) 
(16) 



(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



The terms for the second branc can be obtained anafogously. Using the BPS conditions and keeping only the kinetic 
terms, we get for the Gibbons-Hawking boundary terms 



+ 



Collecting all terms we find 



•^MSA 
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j d^x^/^a^cppsmdcpf 
^ J d^x^ia\cT)UB{a){daf 



/(</.,a)7e(^) + |a^(</.)M^(a<A)^ 



(23) 

(24) 

(25) 



Note that the kinetic term of the field has the wrong sign. This is an artifact of the frame we use here. As wc shall 
show below, it is possible to go to the Einstein frame with a simple conformal transformation, in which the sign in 
front of the kinetic term is correct for both fields. 

The effective action is valid for any model based on STipcrgravity inspired five dimensional models with two brancs. 
In particular it is worth noticing that there is no potential for the moduli. The moduli describe flat directions reflecting 
the no-force condition on boundary branes for BPS-like systems. 

The effective action for the two moduli cj) and cr has a nice interpretation in terms of supergravity. This is expected 
as we started from a two-brane system satisfying BPS conditions. The resulting effective action has a supergravity 
structure. Indeed one can write the Einstein-Hilbert term and the kinetic terms of the moduli as 



d x\/—gi 



where 



/(0,a)7^(4) + + iasSs/s^sa^s 



!($ + $), a=^(S + S). 



(26) 



(27) 



In terms of supergravity, the fields $ and S are the scalar parts of two chiral multiplets. One can read off the Kahler 
potential for these two moduli fields as 
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K = -31n/. 

In the Randall-Sundrum case, we retrieve that 

i^ = -31n(e-'=(*+*)-e-'=(^+^)) 
Notice that in that case one can rewrite the Kahler potential as 



K = -3fc($ + $) - 3 ln(l - e-'=(^+^)), 



where 



(28) 



(29) 



(30) 



(31) 



is the radion superfield. Moreover we immediately sec that, in the Randall-Sundrum case, the field $ can be eliminated 
by a Kahler transformation; this shows that one of the two moduli decouples, leaving only the radion as the relevant 
physical field. We will retrieve this result later when we analyse the equations of motion. 
Coming back to the action above, we redefine the fields in the following way: 



4^ = {l- Aka^ 
0-2 = (1 _ Aka^a) 



.2/3 
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with 



2a' 



Aa' 



then, the gravitational sector can be written as 

1 



•^MSA = 



2A;K|(2a2 + 1) 



(32) 
(33) 

(34) 
(35) 



This is a action of the form of a multi-scalar tensor theory, in which one scalar field has the wrong sign in front 
of the kinetic term. Furthermore, in this frame there is a peculiar point where the factor in front of TZ can vanish, 
namely when tj) = a, which corresponds to colliding branes. We will call the frame in which the action has the form 
of the equation above the bulk frame. 

It is useful to go to the Einstein frame. In order to avoid mixed terms like (9^^)(9''cr), we shall define two new 
fields^: 



= QcoshR, 
a = Q sinh R. 

To go to the Einstein frame we perform a conformal transformation: 
Then [49] 

Collecting everything we get the action in the Einstein frame (where we now drop the tilde): 



S'ef 



1 



2kKl{2a' + 1) 



-9 



6 



12a2 [dQf 

1 + 2a2 g2 2a2 + i 



[dRf 



(36) 
(37) 

(38) 
(39) 

(40) 



Clearly, in this frame both fields have the correct sign in front of the kinetic terms. Note that for a — > (i.e. the 
Randall-Sundrum case) the Q-field decouples. In this case, the field R plays the role of the radion, i.e. it measures 
the distance between the branes. Furthermore, we can identify the gravitational constant: 



167rG = 2fc«;^(l + 2a2). 



(41) 



^Do not confuse the Ricci scalar with the new field R. 
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C. Moduli space approximation: the matter sector 



In the following we introduce matter as well as supersymmetry breaking potentials V{Q, R) and W{Q, R) on each 
branes. We begin with the potentials: to first order in the moduli space approximation we get 

J d'xy/^i [a\4>)V{<t>)\ (42) 

with a4((/)) = ^4/(i+2a^), The expression for a potential W on the second brane is similar with a{<j)) is replaced by 
a{(j). In the Einstein frame we have (dropping the tilde from the metric): 

j rf*ar^/^g-««'/(l+2«')(coshi^)4/(l+2«')F(Q,i^) = j d^x^V,^{Q,R), (43) 

where we have defined 

V^s{Q,R) = Q-^"'/(i+2"'ncoshi?)4/(i+2«')y(g,i?). (44) 
The expression for W{Q, R) in the Einstein frame is 



J d4:cv^Q-«"'/(i+2"')(sinhi?)^/(i+2"'%(0,i?) = J d^xV^W,s{Q, R), 



(45) 



where 

Weff(g, R) = Q-^"'/(i+2"')(sinhi?)4/(i+2«')VF(Q, R). (46) 
For matter, the action has the form 

5W=5W(vl/i,5^«) and S^^^ = S^^^^^g^^'^ (47) 

where denotes the induced metric on each branes and the matter fields on each brancs. Note that we do not 
couple the matter fields \E'j to the bulk scalar field, and thus not to the fields Q and R. In going to the Einstein frame 
we get 

S^^^ =S^^\^i,A\Q,R)g^,) and S^^'^ = S^^^^f^, B\Q, R)g^,), (48) 

In this expression we have used the fact that, in going to the Einstein frame, the induced metrics on each brancs 
transform with a different conformal factor, which we denote with A and B. Wc have neglected the derivative terms 
in the moduli fields when considering the coupling to matter on the brane. They lead to higher order operators which 
can be easily incorporated. The energy-momentum tensor in the Einstein frame is defined as 



^(1) ^ , 1 SS^J^\^,AHQ,R)g) 

with an analogous definition for the energy-momentum tensor for matter on the second brane. 
In the Einstein frame, the total action is therefore. 



(49) 



■ 12a^ {dQf 6 , 
l + 2a2 g2 20^ + 1^ ' 



- j d^x^{VMQ,R) + WMQ,R)) + S^r^H^i,A^{Q,R)g^^) + Sl^\^2,B^{Q,R)g^^). (50) 



The theory derived has a form similar to the one discussed in [50,51], although we have here two scalar degrees of 
freedom in general. However, these scalar fields couple differently to the two matter types on each individual branes. 
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D. The field equations in the Einstein frame 



Prom the action derived above we can now derive the field equations for the metric and the two scalar degrees of 
freedom in the Einstein frame. They are 



+ 



6 



1 



□Q 



UR = 8ttG 



2a2 ^ 
l + 2a2 
6 

l + 2a2 
12a2 



OR 



OR 



12a2 1 
1 + 2q!2 Q2 



/2)j.(2) 



dQ dQ 



(51) 
(52) 
(53) 



In these expression T^*) is the trace of the energy-momentum tensor for each brane's matter; the coupling functions 
ttg^ and a^' are defined as 

a« = ^, a^^) = ^; (54) 

(55) 

We will give the expressions for these quantities in the next section and discuss solutions to the field equations in 
a cosmological setting in Section 4. Due to the coupling of the moduli to the branes, cosmological matter is not 
conserved but satisfies 



8' = 


d\nA 


- 






dQ ' 


dQ 


S." = 


d\nA 




d\nB 




OR ' 


OR 



a}^{d''Q)T, + a}jid"R)T, 



(56) 



for each type of matter i = 1, 2. In deriving this equation we have assumed that the matter fields obey the equation 
of motion. Notice that matter non-conservation is directly linked to the space-time variations of the moduli fields. 

III. COUPLINGS TO MATTER AND STRONG FIELD CONSTRAINTS 

In this section we will be interested in the coupling of matter living on both branes to the bulk. In particular we will 
pay particular attention to the constraints imposed by strong field limits. This, in turn, leads to stringent constraints 
on the parameter a and the allowed values of the field R. 

Let us first notice that matter on the branes couple to the induced metric. The form of the induced metric implies 
that we have to deal with a bi-metric theory. This is a general result for moving branes. In the following we will only 
be interested in the non-derivative couplings between the two moduli and the matter fields. Indeed we are dealing 
with a low energy expansion and derivative interaction lead to higher order operators. If need be these derivative 
interactions can easily incorporated in the following. 

First of all couplings to the two branes results from a Lagrangian 



J ^-9BCm{4>m,9B), 



(57) 



where can either be a scalar, fermion or vector boson field. We will study each case separately. Let us start with 
a scalar field, i.e. the Lagrangian is given by 



= ■^9'^d^(j)Bdv(j)B - V{(t)B)- 

We first write the boundary action in terms of the Einstein frame metric g = Cl~^g leading to 

^g{a''{cj>)n''r''d^,cj)Bd,cj)B - a'{(t>)n''V{(t>B)) , 



I 



(58) 



(59) 
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where fi^ = 1/ f{(j), a). Now we can redefine the scalar field (j) = a{(j))Q.(j)B in such a way that the action becomes 



/ 



-~g \r''d^ct>d.<^ - a''{4>)n^V ( ) ) , (60) 



up to derivative interactions. Let us apply this to a simple renormalizable potential 

V=^m'cPl+Xcl>\ (61) 

In the Einstein frame, the couplings become 

m = a{(f))Q,m, A = A. (62) 

Clearly, in this frame the mass of the scalar field becomes time-dependent. 
We can apply the same technique to vector bosons with action 



5'm = -^ / '^~9b9b'9'bPwP'^> 



(63) 



We find that in the Einstein frame 

9E = 9, (64) 

as expected vector bosons are conformally coupled and therefore the gauge coupling constant are time- independent. 
Finally let us deal with Dirac fermions: 

•Sm = y V-5b(#7b-C)^'0 - m-0V'), (65) 
where 7^ are the gamma matrices for the induced metric. Using the results derived in [52] , we can rewrite the action 



as 



^m = J (i^rD^^i' - ma{4,)ni>ii) , (66) 

where we have used the conformality of the coupling of massless Weyl fermions and defined tjj = {a{cj))Q)^^^ip. Again 
we find that fermion masses become time-dependent 

m = a{(j))Q.m, (67) 

Notice that the time dependence of the fermion and scalar masses is the same. This is not true when comparing 

different particle species on the positive tension brane with species on the negative tension branc. However, particles 
on the negative tension brane are candidates for dark matter (they couple to the standard model particles only via 
gravity). Constraints imposed on dark matter by Equivalence Principle violations are less restrictive. 

In conclnsion wc see that, in the Einstein frame, the only time-dependent coiiplings in the matter Lagrangian are 
the masses. Indeed, in a general frame we have the following invariant relating Newton's constant to masses 

I = Gm"^ = Gm". (68) 

which is a dimcnsion-lcss quantity. In particular we find that in the frame where masses are time-independent the 
Newton constant becomes time-dependent. As such / gives the true measure of time variations in the brane- world 
models that we consider. 

The fact that the coupling constants are frame-independent (see eq. (64)) and therefore space-time independent 
means, that in the theory we consider we would not expect a time-varying fine-structure constant, as reported in 
[40,41]. In order to explain a time-varying coupling constant, one needs to couple vector bosons directly to the bulk 
scalar field and thus to Q and R. We do not consider this possibility in this paper. 

Let us now consider the constraints imposed by the strong-field limits. For that it is convenient to write the moduli 
Lagrangian in the form of a non-linear sigma model with kinetic terms 

HM^^^ (69) 
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where i = 1,2 labels the moduli 4> = Q and ^ = R. The sigma model couplings are here 

12^2 16 , , 

Notice the potential danger of the a ^ limit, the RS model, where the coupling to Q becomes very small. In an 
ordinary Brans-Dicke theory with a single field, this would correspond to a vanishing Brans-Dicke parameter which 
is ruled out experimentally. Here we will see that the coupling to matter is such that this is not the case. Indeed we 
can write the action expressing the coupling to ordinary matter on our brane as 

Sm = Smicpm, A^g), A = a(</.)fi, (71) 

where we have neglected the derivative interaction and expressed^ qb = Ag. Notice that A enters in both the coupling 
of matter to the brane and the time variation of masses. As such it represents the actual time-dependence of couplings 
in our models. Let us introduce the parameters 



aQ = dQ In A, or = Or In A. (72) 



We find that (A = 4/(1 + 2q2)) 



leading to 



A 



^ = (coshi?)4, (73) 



a^A 1 Atanhi? 

Observations constrain the parameter 

9 = -I'^aidj (75) 
to be less than 10~^ [53]. We obtain therefore a bound on 

4 tanh^ R 

3 1 + 2a2 + 6(l + 2a2)" ^ ' 

The bound implies that 

a < 10-^ R < 0.2 (77) 

This already rules out the supergravity model based on = 1/12 [44]. The smallness of a indicates a strongly 
warped bulk geometry such as an Anti-de Sitter spacetime. In the case a = 0, we can easily interpret the bound on 
R. Indeed in that case 

tanhi? = e-'=(''-'^), (78) 

i.e. this is nothing but the exponential of the radion field measuring the distance between the branes. We retrieve 
the well-known fact that gravity experiments require the branes to be sufficiently far apart. When a ^ but small, 
one way of obtaining a small value of R is for the hidden brane to become close from the would-be singularity where 
a{(j) = 0. In the following we will analyse the cosmological behaviour of the present model and in particular the 
robustness of the condition i? <C 1 to cosmological evolution. 



^The parameter B which defines the coupling of matter to the second brane is similarly defined hy B = a{a)Q. 
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IV. COSMOLOGICAL EVOLUTION OF THE TWO BRANE SYSTEM 



As we have seen from the last section, the field R has to be small today in order for the theory to be in agreement 
with observations. This required smallness of R could result from one of the following possible features of our model: 
firstly, the minimum of VeS or WeS could be at small values of R, and consequently R would have been driven towards 
this minimum in the very early universe; the form of the potential has to be derived from an underlying theory, whose 
form is unknown. A much more interesting alternative is that R could be driven towards small values during the 
late cosmological evolution, i.e. after nucleosynthesis, via an attractor mechanism. In scalar-tensor theories such a 
mechanism is well known (see e.g. [37-39]) and the question we address in this section is whether such an attractor 
mechanism exists for the brane world models we have discussed so far. Clearly, in our model there are two moduli 
fields, Q and R and even if R is rapidly driven towards small values, there is the danger that the dynamics of the Q 
field has no attractor behaviour. This could jeopardize the cosmology of the model we consider. Furthermore, even 
if such an attractor mechanism exists, it is not a priori clear whether it is efficient enough. 

To address these important issues, we now discuss cosmological solutions to equations (51), (52) and (53). In this 
case the field equations are as follows (a 7^ 0): The Friedmann equation reads 

H'=^-^ iPi +P2 + v.. + W.,) + + ^R^ (79) 

where we have defined Q = exp cj). The second Einstein equation is 

H + H^ = (pi + 3pi +P2 + 3p2 - 2Vef{ - 2Weff) 



4a2 



l + 2a2^ 1 + 2q!2 
The field equations for R and ^ read 



^2 if (80) 



R + 3HR=-&TrG^-^^^ 



OR ' dR (^i~^^'i)+"V(P2-3p2) 



(81) 



12q!2 

For the model we consider here, we have 



dVeS , dWes (1), , , 

-^T- + -^J- + "rf. (Pi - 3pi) + < (P2 - 3p2) 



2a2 



a 



(1) 



1 + 2a2 ' 
tanhi? 



ai^) = - 



2a^ 



1 + 2a2 ' 



l + 2a2 



(2) _ (tanhi?)-! 



l + 2a2 



(82) 

(83) 
(84) 



We will study numerical solutions of the system in more detail below, but if we consider matter on the positive 
tension brane only we can draw some conclusions concerning the evolution of the fields 4> and R: 

We considerer that throughout the radiation- and matter-dominated eras the potentials Vgg and Wefj are negligible. 
This is a similar condition as that imposed in quintessence models. Then, in the radiation-dominated era, the traces 
of the energy-momentum tensors vanish, implying that 



R- 



(85) 



in such a way that R remains small if it is small initially. Moreover, there is no change to the time dependence of the 

scale factor. 

In the matter-dominated era, the equations of motion read 



R + 3HR 



8nG 



N 



Pl + 



2a' 



12 



inG 



N 



PiR' 



3(l + 2a2) 



+ 



1 + 2q!2 



-R, 



I + 3H<p = — pi 



2 
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together with the conservation equation 



pi + ZHpi « -2q!>i(/' (86) 

where we assume that both a and Rq are small to comply with the tests of the Equivalence Principle. Assuming that 
the energy density of matter dominates the expansion rate, the solutions to these equation is given by 

/ ^ \ -3-2aV3 / ^ X 2/3-4aV27 



together with 



^ = ^o + 3ln-,i? = i^.^-J +R.^-] (88) 



as soon as t » to- Note, that this solution is consistent with our assumptions on the expansion rate as long as R is 
small at t = to- Notice further that R converges towards zero, hence retrieving general relativity at late times and 
that the change of the expansion rate is of order and therefore small. The solution above is in good agreement with 
the numerical solution described below. We therefore conclude that during the matter-dominated era R is driven to 
zero, provided R was small at the beginning of the matter era. 



A. Numerical Results 



We now describe numerical solutions of the cosmological equations. There are different interesting cases to study. 
Firstly, all matter types could live on the positive tension brane; a more intriguing alternative might be that the dark 
matter lives on the negative tension brane, whereas radiation and baryons live on the positive tension brane. We will 
discuss these cases below. 



1. No Potentials: 



We begin the analysis with the assumption that the potentials V and W are identically zero (in other words, that 
the brane potentials are unbroken from their supersymmctric values); thus, we follow the evolution of the fields during 
radiation and matter domination. Both matter and radiation live on the positive tension brane. The calculations are 
made with a = 0.01. 



1 

0.8 
0.6 
0.4 
0.2 
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FIG. 1. The evolution of R with different initial conditions for the case of radiation and matter on the positive tension brane 
and no matter on the negative tension brane. We find that R is driven towards zero. 



The evolution of R and (p are shown in fig. 1 and fig. 2, respectively. One can clearly see that during radiation 

domination both fields are frozen in, because the trace of the energy momentum tensor is effectively zero. Soon after 
matter becomes important, both fields are forced to evolve due to the non-vanishing trace of the matter energy- 
momentum tensor. For the initial conditions we have chosen the constraint R < 0.2 can be fulfilled. However, if R is 



11 



phi 
2 

1.5 

1 

0.5 



10^ 10^ 1 



FIG. 2. The evolution of (ji for different initial conditions with the same cosmology as in Figure 1. 

initially large, it has not enough time during matter domination to evolve to small enough values. Thus, there is a 
constraint that R has to be small enough during radiation domination. 

If wc put matter on the second branc as well, the evolution of R is modified, as can be seen in fig. 3 (sec fig. 4 
for the evolution of 0). For these cases, the matter on the second brane does never dominate the expansion of the 

universe. However, due to the coupling function a)j , the evolution of R and (j) are affected strongly by matter on the 

second branc. Indeed, as can clearly be seen, R is driven faster towards small values even with only a small amount 
of matter on the second brane. The field (j) stays longer constant after radiation domination if there is matter on the 
second brane. 



R 




1 +z 



FIG. 3. The evolution of R with different initial conditions for the case of radiation and matter on the positive tension brane 

and pressureless matter on the negative tension brane. We find that R is driven towards zero for all cases. Note, that if the 
ratio P2/P1 grows, R evolves faster towards zero. The ratios of P2/P1 are given by 0, 0.11, 0.25, 0.42 and 0.5, for the curves 
from the right to the left. 

We conclude that there is an cosmological attractor which drives R towards small values and that the efficiency is 
enhanced if the energy density of matter on the second brane is not negligible initially. 



2. Including Potentials: 

We now include potentials on the branes. We will study the case where there is only a potential on the positive 
tension brane. For the potential we do assume that it starts to dominate the energy density only recently in the 
cosmic history, because we are interested in the evolution of the moduli fields R and (j). In principle, both and R 
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FIG. 4. The evolution of (ji for different initial conditions with the same cosmology as in Figure 3. 



arc candidates for the dark energy. However, we do not address the coincidence problem, nor do we want to provide 
a model for dark energy in this paper (for a discussion about dark energy arising in particle physics models, see e.g. 
[56-64] and references therein). Instead, we use the simplest model for the potential, namely we assume that V has 
the form of an exponential potential. To be specific, we assume that supersymmetry is broken by tuning Ub away 
from the form (8) by setting 

V={T- l)Ake°"^. (89) 

Here, T 7^ 1 is a supersymmetry breaking parameter (1/) is the bulk scalar field). Expressed in terms of (j) and R we 
have 

F((/., R) = 4(r - i)A;e-i2«'*/(i+2"') (coshi?)(*-^"')/(i+'"') . (90) 

Notice that for R close to zero, this is nothing but an exponential model with the field (p playing the role of a 
quintessence field. In the following we will set Ak{T — 1) such that the universe starts to accelerate at a redshift 
around 1. 

Omega 




FIG. 5. Evolution of the density parameter fij as a function of redshift for radiation, matter and the scalar fields. When the 
potential of the moduli dominates, the universe is accelerating. Note, that in order to explain the values for the energy density 
of dark energy, one has to fine-tune the parameters of the theory. For these plots we have set a = 0.01. The dark matter lives 
on the positive tension brane, there is no matter on the negative tension brane. 

The evolution of the cosmological parameters is shown in fig. 5. After the usual matter dominated era, the universe 
becomes dominated by the potential energy of the fields and starts to accelerate. 

The evolution for R is shown in fig. 6. We have chosen the same initial conditions for R as in fig. 1 and it can be 
seen that the evolution of R is not much affected by the presence of a potential which dominates today. In fig. 7 we 
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FIG. 6. Evolution of the field R as a function of redshift with different initial conditions. The constraint today is R < 0.2. 

have plotted the evolution of cj). Note that as soon as the potential dominates the expansion, the evolution of (j) is 
affected. 

phi 
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FIG. 7. Evolution of the field as a function of redshift with different initial conditions. Note that when the potential energy 
starts to dominate, 4> changes its time-evolution. 



3. Time-variation of Masses or the Gravitational Constant 

As explained in Section 3, the quantity /, defined in eq. (68), is an invariant under conformal transformations. 
This quantity specifies the variation of masses in the Einstein frame (or the gravitational constant in the frame where 
the masses on the positive tension brane are time-independent). The variation of the gravitational constant (and 
therefore of J) is constrained by nucleosynthesis (see e.g. [54] and [55] for a recent discussion). We have plotted the 
evolution of / for different initial conditions, one with initial conditions leading to an evolution in agreement with 
observations (fig. 8) and one extreme case which is ruled out by nucleosynthesis constraints (fig. 9). Therefore, 
the initial conditions for the fields R and 4> are not arbitrary, but constrained by nucleosynthesis. The details of the 
evolution of / is strongly dependent on the matter contents on the branes. 

B. The Five— Dimensional Picture 

Having applied the four-dimensional effective theory to cosmology, we turn now our attention to the five-dimensional 
interpretation of the solutions above. Specifically, we consider the brane positions, which are related with (p and R 
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FIG. 8. Evolution of //Jo (eq. (68)) for initial conditions which are allowed by nucleosynthesis. The value of the gravitational 
constant at nucleosynthesis was of order one percent larger than its value today. 
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FIG. 9. Evolution of I /lo (eq. (68)) for initial conditions which are forbidden by nucleosynthesis. 

via eqns. (32), (33), (36) and (37). They are shown in fig. 10. It becomes clear from eq. (37), that R = Q corresponds 

to cr = l/Aka^, i.e. the negative tension brane is attracted towards the singularity. Negative vahies of R do not make 
sense in this description, when the transformations (32), (33), (36) and (37) are applied. Therefore, cosmological 
solutions based on the moduli space approximation in which R is negative do not have a sensible five-dimensional 
interpretation. 

There are two obvious interpretations of this result: 

• At the singularity (i.e. at i? = 0) the negative tension brane will be destroyed. There is only one remaining 
scalar degree of freedom ^. The late time evolution is essentially that of a one brane system. 

• The negative tension brane is repelled from the singularity. It is tempting to speculate that the repelling can 
be described by an effective potential in the effective four-dimensional description. This option will be briefly 
discussed below. 

Thus, for a complete understanding of brane cosmology with a bulk scalar field and two branes one needs a better 
understanding of the singularity, which might appear from string theory (for a discussion on singularities in string 
theory sec e.g. [65] and [66]). Nevertheless, we would like to emphasize that the net effect remains the same: R is 
initially attracted towards small values. 
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FIG. 10. Evolution of the brane positions for the cosmology as in figure 5. 

C. Avoiding the Singularity? 

Given that during the matter dominated epoch the negative tension brane is attracted towards a singularity, we 
will now investigate the response of the system with the adoption of the second option mentioned above; i.e. we add 

a potential which drives the brancs away from the singularity. What arc the conditions for such a potential? First, 
whatever type of matter wc have on the brancs, near the singular point, i.e. Q — {(f> ~ — oo) and/or i? = 0, the 
potential term in the Klein-Gordon equations must dominate. This can be achieved by an inverse-power law potential. 
A second ingredient is that for consistency with the limit a ^ 0, the derivative of the potential must be proportional 
to or some higher power (see eq. (53)). A simple example of a potential which fulfills these requirements is 

Kdd = Vo{Q''"r)-\ (91) 

However, once the negative tension brane is repelled from the singularity, the field R grows again and is at the 
present epoch, generically, too large to be consistent with observations. Therefore, adding a potential which repels 
the branes from the singularity jeopardizes the valuable properties of the attractor mechanism during the matter 
dominated epoch. This can be easily understood: the attractor mechanism is due to the coupling functions and 

Q!j^^ appearing in the Klein-Gordon equation (81). If at some point during the matter dominated epoch these terms 
become negligible compared to the potential term, the field R will be driven by the gradient of the potential. For the 
potential above, this will be towards larger values of R. Then, at some point the potential term in the Klein Gordon 
equation becomes negligible again compared to the matter-coupling terms, which drive R towards smaller values until 
the potential term dominates again. During this interplay between the terms in the Klein-Gordon equation, R will 
be driven towards larger and larger values because the density of matter drops and it takes more and more time for 
the potential term to be negligible. The expansion of the universe is never dominated by the potential energy of the 
moduli fields for the potential (91). 

What happens in the dark energy dominated regime? This depends clearly on the nature of the dark energy. If the 
fields (p and R are responsible for dark energy, then we need to modify the potential at large R, as otherwise R will 
be too large to be consistent with observations and, as explained above, the potential (91), for example, will not lead 
to a accelerating universe. An example which delivers an accelerating universe would be the potential (90) added to 
(91). However, this corresponds to the situation where the total potential has a minimum at -C 1, which has to 
be fine tuned. A cyclic universe can be modeled by the inclusion of a potential which is negative in value for some 
range of the moduli fields [48] . 

If instead the dark energy has a different origin and is a field which is confined on a brane, then there is an important 
interplay between the potential term and the coupling term in (81). We will investigate this possibility in future work. 

V. CONCLUSIONS 

In this paper we have investigated the cosmological evolution of brane world moduli for a general class of brane 
world models with two boundary branes and a bulk scalar field. The theory we have used contains the Randall- 
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Sundnim model as a special case. The parameter a regulates how strongly the bulk geometry is warped; small values 
of a correspond to highly warped geometries, whereas large values of a correspond to only slightly warped geometries. 

We have obtained the four-dimensional effective theory from the five-dimensional theory using the moduli space 
approximation. As we are interested in late time cosmology, namely from nucleosynthesis on, the moduli space 
approximation should be accurate enough. We have discussed observational constraints and found that the parameter 
a must be small, pointing towards a warped bulk geometry. Furthermore, the field R must be small enough today, 
meaning that the distances between the branes must be large enough. 

As long as matter is not directly coupled to the bulk scalar field, and therefore to (p and R, the theory predicts that 
coupling constants such as the fine structure constant are not space time dependent. Thus, if the observations made 
in [40] and [41] hold, more complicated models of brane worlds have to be investigated. 

Finally, we have discussed the cosmological evolution of the moduli fields. We focused in particular on the question of 
whether R is driven towards small values during the cosmological evolution in order to be consistent with observations 
today. Our findings indicate that there are cosmological attractors for the field R, similar to the attractor found in 
scalar-tensor theories [37-39] . We have also found that the efficiency of the attractor strongly depends on the matter 
content on the negative tension brane. In the five-dimensional picture, the attractor solution drives the negative 
tension brane towards a bulk singularity. When R = 0, the negative tension brane hits the singularity. Negative 
values of R correspond to the situation where the singularity is between the branes. As the four-dimensional solution 
suggests an oscillating behaviour of R around 0, one has to be careful when interpreting this behaviour. Indeed, a 
proper understanding of the singularity is absolutely necessary in order to fully understand the cosmology of the two 
brane system. In the case with the brane repelled by the singularity we have argued that in this case the details of 
the evolution of the two brane system can be significantly altered. 

We would like to emphasize again that the class of brane world models we have considered in this paper predict 
the coupling functions (83) and (84). These functions control the evolution of the fields <!) and R and the coupling to 
the different matter species. 
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